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Abstract 

CZ2 . 

O ' We draw attention on the fast oscillatory deviations from the residual 

" 5j ' non-ambipolarity in the case where electrons are driven by strong mag- 

netic perturbations. 

' 

O ■ If the turbulent diffusion of the electrons and ions are different one can 

'— 1 \ expect that a radial electric field develops (such that the ambipolarity is rein- 

| ■ stored) and that poloidal plasma spin-up can occur. However it has been shown 

that the diffusion of particles in turbulent fields is ambipolar. In Ref. |lj] the dif- 
X/"} • ference in the radial fluxes generated by fluctuating fields has been considered 

\ and the role of the ion polarization drift has been underlined. The radial ion 

■ flux partially cancels the radial electron flux and the ambipolarity is reestab- 
lished approximately, to the order l/e±, which is a very small quantity (e± is 

■ the plasma transversal permitivity constant). The problem has been examined 
in detail for instabilities having a substantial magnetic component in Ref. 

""^ ■ When the unstable modes that drive the radial transport are bounded inside 

O ! plasma (i.e. it can be assumed that there is no exchange of momentum with the 

* ^5 ' region outside plasma) the fluxes are ambipolar in average over a spatial region 

including several mode-resonance surfaces. Pointwise non-zero radial currents, 

■ even reduced by the fraction l/e_i_, are able to build up charge filaments and 
can cause the plasma to be Kelvin Hclmholtz unstable. 

We return to this problem to examine a situation where the non-ambipolarity 
is switched on by an external source acting on the electron component. We 
assume that the electrons are radially driven by a magnetic perturbation which 
imposes a constant radial flux starting at t = 0. The ion polarization drift 
re-establishes the equality of the fluxes (to order l/e_i_) but in this particular 
case a stationary equilibrium cannot be reached. The existence of a radial 
non-zero current yields a non-zero time derivative of the radial electric field, 
i.e. a non-zero time derivative of the poloidal velocity. The plasma would be 
accelerated without limit in the poloidal direction, but the torque competes with 
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the magnetic pumping, a very effective damping mechanism. We draw attention 
to the oscillatory variation of the velocity, which occurs an a scale much faster 
than the dissipative decay. 

Such "events", consisting of sudden creation of non-ambipolar fluxes, fol- 
lowed by a fast plasma response in view of reinstoring the ambipolarity, can 
appear in a random space and time sequence and on the average can affect 
the plasma dynamics. To look in detail to only one event we shall adopt an 
"initial value" point of view. We consider the slab geometry with x the radial 
coordinate increasing from the reference point toward the centre of plasma, y 
the poloidal coordinate, z is directed along the shearless magnetic field such as 
the mixed product of the three versors is positive. The radial electric current is 
switched on at t = , j x (t = 0) = —env% (we note e = |e|) and for simplicity we 
shall assume that v e x is constant in time and uniform in space during this single 
event. From the momentum conservation we find (with c and da respectively 
the light speed and the Alfven speed), 



We have neglected the damping and the diamagnetic flow. From the initial 
non-ambipolar current {—env% — j x (0) ^ 0) it results a time-growing poloidal 
velocity whose magnitude is inverse proportional to the very large factor repre- 
senting the plasma transversal permitivity e± = 1 + (? jv\. 

Since the time-derivative of the poloidal velocity is constant, the plasma 
rotates indefinitely in the poloidal direction with higher and higher velocity. 
Naturally this requires the consideration of saturation mechanisms and makes 
desirable a time-depending investigation. We write in more detail the ion mo- 
mentum equations on x and y 




(1) 




(2) 




(3) 




(4) 



The charge conservation gives 



= j x + £q 
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where 



jx (t) = -env% + env l x 



(6) 



2 



From these relations and using Eqs. 
for the poloidal ion velocity. 



and (pNq) we obtain a single equation 
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where £q is the vacuum permitivity, (i) is the Heaviside function and 



H 
dx 



This equation contains higher time and space (x) derivatives. In order to 
simplify the problem we assume that the externally imposed electron flux is 
uniform on the radial direction (for the region of interest) and this renders the 
entire model invariant to translation in x coordinate. Then we shall restrict 
this discussion to solutions which are independent of x. We shall neglect the 
problem of the stability of the uniform solutions at perturbations in x. Then 
we have 

-v e x @{t) = (b 1 +b 2 )w + b 3 w 

Here 
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The quantities appearing in the equation have been non-dimensionalized: t — > 
t' = t/r, y — > y' = y/l, v — > v' = v/vq, where we take as units: r = £1^, 
I = p s the ion Larmor radius at the electron temperature, vq = (3v^ h where [3 
is the radial magnetic perturbation B x normalised to the main magnetic field, 
= B x / B. v^ h is the electron thermal velocity. Finally we remove the primes. 
The dots means derivative to the time variable. Then the coefficients are 
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The solution is obtained by Laplace transform and it reads 
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The initial conditions are given for the first and the second derivatives of the 
velocity, respectively wq, and w\. In physical units 
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The frequency of the oscillations is in physical units 
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Now we include the effect of the damping due to the magnetic pumping, 
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where v is the appropriate decay rate 0, ||, B. The 



equation restricted to the time domain is now 
-v%Q{t)=d Q vi + d 1 (^j + d 2 ^ 

where 6 (t) is the Heaviside function and 
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The solution of the equation is obtained by the Laplace transform and is written 

i \ ^ / \ d-i v\c , . 

u ft) = w B m + fig-r" "! ex P (7*) + 

do 7 
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exp (a r t) [(a r a r + ajO;,) cos (a^i) — (a^ov — a r oti) sin (a^t)] 



Here a, /3 = a* and 7 are respectively the two complex conjugate and the real 
roots of the polynomial equation s 3 + {d 2 /d 3 ) s 2 + (di/d 3 ) s + (do/d 3 ) — 0. We 
have for the real (r) and imaginary (i) parts: 
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The notations are 
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and a = 1 — b 



(I2V1 + d 3 v 2 



d 3 v 1 



d v 



d 3 vi 

The initial order i-derivatives of the velocity are noted Vi, i = 0, 1, 2. They 
are taken to reflect the sudden onset of the electron flux, i.e. the initial poloidal 
velocity is zero, vq = v y (t = 0) = and the first derivative is obtained from the 
assumption that at t = the ions did not yet moved radially (v l x 
means that 



0) which 



at 



<) 



The second derivative is taken zero, v 2 = 0. The asymptotic value of the poloidal 
velocity is insensitive to this parameter. 



1 Discussion 

The most interesting aspect of this solution is the oscillatory behaviour which 
is strongly controlled by the non-dissipative part of the equation (i.e. terms not 
depending on v). The dynamic is very clear: just at t — the ions acquire a 
high time derivative of their poloidal velocity, imposed by j x (t = 0) = —env%. 
The velocity v l y starts from zero and grows very fast to values which can be 
higher than those required to cancel the non-ambipolar current j x . Then they 
reverse the current and since the electrons are constrained to follow the external 
drive (magnetic radial perturbation), the tendency of rotation reverses and so 
on. The decay due to the magnetic pumping appears on longer time scales and 
this indeed reduces the plasma rotation, without completely eliminating the 
oscillations. A plot, Fig. ([!]) of this evolution is shown, on very large time scale 
(useful to see the atenuation). 

The asymptotic value of the poloidal velocity is (with our choices of initial 
conditions) 

4 {t^l- l ) = {-v%)enB z v- 1 

We can see that the nature of these oscillations is similar to that caused by a 
deviation form electric neutrality, when the frequency is the plasma frequency 

1/2 

uj p = (e 2 n/eom e ) . The important effect of the of these oscillation is in the 
energy balance. The plasma rotation takes energy from the source of the initial 
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Figure 1: Time variation of the poloidal velocity 



magnetic perturbation and from the electron kinetic energy, since both support 
the radial electron current. This energy is dissipated via the collisional ion 
viscosity represented by the magnetic dumping. A more complex treatment of 
this problem should include a detailed electron dynamics, under the effect of 
the external magnetic perturbation (which must be quantified). However, the 
presence of an oscillatory behaviour is still expected and suggests to take this 
into account when the efficiency of magnetic stochasticity-induced transport is 
considered. 
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